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Abstract 

Inflation of cosmic gauge and global strings is investigated by numerically 
solving the combined Einstein and field equations. Above some critical 
symmetry-breaking scales, the strings undergo inflation along the radial di- 
rection as well as the axial direction at the core. The nonsingular nature of 
the spacetimes around supercritical gauge and global strings is discussed and 
contrasted to the singular static solutions that have been discussed in the 
literature. 
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I. INTRODUCTION 



Cosmic strings are linelike topological defects that may form as a result of a phase 
transition in the early universe. If a string is associated with a magnetic field, it is called a 
gauge string, otherwise it is a global string. They have attracted much attention because of 
their cosmological importance: deficit angle in the spacetime geometry and a candidate for 
the seed of structure formation in the early universe. 

It was proposed that topological defects can inflate if the symmetry-breaking scale sat- 
isfies T] ^ T] c ~ 0(m p ) in Refs. [jM. This was later verified in numerical simulations by 
Sakai et al. ||. They found, in particular, that the critical value of i] for domain walls and 
global monopoles is r] c ~ 0.33m p . Then what about cosmic strings? There is no reason 
why we exclude cosmic strings out of the topological inflationary category. Recently, it was 
numerically proved that a (2+1) dimensional gauge string can inflate by de Laix et al. [|J. 

In this paper, we shall numerically solve the combined Einstein and field equations for 
a gauge and a global string in (3+1) spacetime dimensions. For the gauge string, we find 
that the core inflates if rj > 0.25m p with unit winding number in the critical coupling case 
(Bogomol'nyi limit). For the global string, 7] > 0.23m p . The critical values decrease as the 
winding number increases. For the gauge string, the critical value also decreases slightly as 
the coupling of the gauge field to the scalar field becomes weaker than the self coupling of 
the scalar field. 

The asymptotic spacetime of a gauge string is known to be conical [f| . This spacetime 
exhibits a deficit angle A = 8ttG[i, where /i ~ rf is the mass per unit length of the string. 
When the symmetry-breaking scale is sufficiently large, the deficit angle exceeds 2tt and 
analyses of the static solution show that the spacetime possesses a physical singularity 
outside the core of the string However, from our numerical results we know that 

supermassive strings are dynamical and undergo inflation at the core. Therefore, we believe 
that the static treatment of supermassive strings loses its validity and that we should treat 
them in a time- dependent way. 
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For global strings, the singularity exists regardless of the symmetry-breaking scale. Many 
people have tried to find a static solution of a global string and they found that there also 
exists a physical singularity outside the core of the string P-fT^]. What was suggested to 



remove this singularity is again a time-dependent treatment of the string. Gregory (13 
introduced a specific metric which has an axial time dependence and showed that this 
spacetime is nonsingular. 

In our work, we follow the evolution of supermassive gauge and global strings in a general 
time-dependent metric and show that no singularity develops in the spacetimes around the 
strings. 

In the next section, we solve the Abelian Higgs model of a gauge string and discuss its 
inflation and spacetime geometry. Sec. Ill is devoted to global strings. Our conclusions are 
summarized in Sec. IV. In Appendix, we show the equations in detail and the numerical 
algorithms. 



II. GAUGE STRING 

Let us consider the Abelian Higgs model of a string with a complex scalar field $ and a 
£7(1) gauge field A^, coupled to gravity. The action is 

where 71 is Ricci scalar and the Lagrangian C is 

C = -D^D^ - ^F^F^ - V($<5) . (2.2) 

Here the covariant derivative is = V M + ieA^, the gauge field strength is F^ = V fl A u — 
VvA^, and the potential is 

where e is the coupling constant between gauge and scalar fields, A is the self-coupling 
constant of the scalar field, and r\ is the symmetry-breaking scale of the scalar field. We use 
a general cylindrically symmetric metric 
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ds 2 = -dt 2 + B(t, r) 2 dr 2 + C(t, r) 2 r 2 dQ 2 + H(t, r) 2 dz 2 . (2.3) 

(We use ti — c — 1, G — 1/m 2 .) In these coordinates, we can use a time-dependent 
generalization of the usual Nielsen- Olesen anzatz for the scalar and gauge fields 

$ = 0(t,r)e in< \ 



n 



= a(t,r)V^, 

e 

where n is the winding number of the string. We also set the usual boundary conditions 

<P(t,0) = 0, 0(t,oo)=77, 

a(*,0) = 0, a(t,oo) = l. (2.4) 
For the above action, the Einstein's equation is 

G/mu = 87rGTfj, u , 
where the energy-momentum tensor is given by 

T^ v = D^D V § + D 1/ <$>D t j& + g af3 F^ a F u(3 + g^C . (2.5) 
The scalar field equation is 

dV(<5>$) 
D = - - 

and the gauge field equation is 

The detailed equations are given in the Appendix. 

The static solution at a large distance from the core is 

d s 2 = -dT 2 + dR 2 + (1 - 4Gfx) 2 R 2 d6 2 + dZ 2 , (2.6) 

where \i is the mass per unit length of the string and proportional to rf. The metric ( |2.6| ) 
exhibits a deficit angle A = 8irGfi. As rj increases, the deficit angle also increases. When 



the deficit angle exceeds 2%, the static solution ( |2.6|) ceases to exist and we expect that the 
string becomes dynamical and undergoes inflation at the core in this regime. 

Let us consider the role of the winding number for inflation. The coupling between scalar 
and gauge fields is given by 

m q \ 2 A 



P -{mJ 2e 2 ' 

where m s and m v are masses of scalar and gauge fields. In the limit of critical coupling 
(Bogomol'nyi limit), (3 = 1, the mass per unit length is fi = 2ir\n\r] 2 |14|]- Then the critical 



value of 7] at which the static solution (|2.6| ) ceases to exist is proportional to 1/ \fn. Therefore, 
as the winding number increases, the core of the string can inflate at lower symmetry- 
breaking scales. Another way to explain the inflation of topological defects is the core 
size being bigger than the horizon scale Under this condition the gravitational effect 

becomes important in the core region and the defects are supposed to inflate. The energy 
density of the string can be roughly estimated by 



2 P _ n V 



6 2 6 2 ' 

where 6 is the core size and H is the Hubble parameter. The core size is then 

j^~V^V- (2-7) 

The core of the string can inflate when its size becomes comparable to the horizon scale, 
~ 1. This condition gives again r] c ~ horn Eq. ( [2.7|) . In other words, for a given 

symmetry-breaking scale rj, the core size S is proportional to ^fn. Therefore, as the winding 
number increases, the wider core of the string can inflate at lower symmetry-breaking scales. 
Figure [I] shows the core sizes as functions of the winding number for several 77's in the flat 
spacetime. 

One important effect of the gauge field coupled to the scalar field is that it makes the 
string well localized and, thereby, suppresses the energy divergence along the radial direction. 
As the coupling constant j3 increases, the effect of the gauge field becomes weaker, therefore, 



the core size becomes bigger. The string has a better chance to inflate with the bigger core 
acquired as (3 increases. 

We numerically solve the field equations introduced earlier. Initially we assume a flat 
spacetime and obtain the initial scalar and gauge fields, 0(0, r) and a(0, r), by numerically 
solving the flat-spacetime field equations. We then solve the time-dependent field equations 
to trace the evolution of these fields as well as the gravitational field g^ v . We define the core 
size by the proper radius Cr at which the scalar field is = rj/2 and determine whether 
inflation occurs by observing the growth of this proper radius. For n = 1 and (3 = 1, the 
critical value of r\ for inflation is found to be rj c ~ 0.25m p . This value is remarkably close to 
the critical value 0.255m p found by Laguna and Garfinkle ||, at which the static solution 
becomes singular. The critical value decreases as n and (3 increase as expected. The critical 
values of 77 are given in Table | for several n's and (3 , s. Figure ^| shows the evolution of the 
scalar field configuration as a function of the proper radius, t and r are scaled in the unit 
of the horizon size at the center, Hq 1 = [8irGV((j) = 0)/3] _1//2 . We can see clearly that the 
proper radius grows very rapidly (exponentially) in the core region. 

Having found the inflationary behavior in the supermassive stringP], it is worth while 
to discuss its spacetime. For the supermassive string, the static solution ( |2.6p is no longer 
valid. However, since the energy-momentum tensor of the string falls off rapidly outside the 



core, the asymptotic form of the metric must be one of the two Levi-Civita metrics ||15|| , 

ds 2 = -dt 2 + dz 2 + dr 2 + ( ai r + a 2 ) 2 d9 2 (2.8) 

which is conical, or 

ds 2 = (hR + 6 2 ) 4/3 (-rfT 2 + dZ 2 ) + dR 2 + (b t R + b 2 y 2/3 d8 2 , (2.9) 



which is a special case of a Kasner metric |fL6|| . The qualities of the metrics (|2.8|) and ( |2.9| 



depend on the signs of the constants a l5 a 2 , b\, and b 2 - When a\ and a 2 have the same 



1 "Supermassive" was named for the gauge string to which the static solution ( |2.6[ ) no longer 
applies 
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sign, the metric ( |2.8|) is equivalent to the metric ( |2.6|) . For b± and 62 with different signs, 
Laguna and Garfinkle || examined that the asymptotic spacetime of the supermassive string 
is described by the metric ( |2.9| ). For a\ and a-i with different signs, Ortiz |7|] demonstrated 
that the metric ( |2.8| ) can also be a solution. In both cases, the spacetimes possess physical 
singularities at r = — a 2 /ai and R = — &2/&1- The singularity approaches the core as 77 
increases. Nevertheless, it is not clear whether it is physically sensible to have a singularity in 
the string spacetime. Since we have shown that the supermassive string undergoes inflation, 
we believe that the static treatment is no longer valid and the string should be treated in a 
time-dependent way to describe its dynamical nature. We expect that the time-dependent 
treatment removes the singularity from the spacetime of the supermassive string. 

One simple way to examine the existence of a physical singularity is to analyze the 
metric. The Riemann tensor contains the second derivatives of the metric terms. If the 
metric terms do not contain cusps and are smooth functions of the coordinates up to the 
second derivatives, then the scalar invariants calculated by the Riemann tensor are finite 
and it is safe to say that no singularity develops in the spacetime. While we perform the 
numerical calculation, we follow the evolution of the metric. Figure [5] shows logarithmic 
values of the metric terms B, C, and H at several moments of time. They are smooth 
functions of r and change also smoothly in t. We also calculate the Kretschmann scalar 
RaPisR a ^ & '■ This invariant is finite everywhere and every moment in time as shown in 
Fig. |j. Therefore, we conclude that the time-dependent treatment of supermassive strings 
makes a singularity-free spacetime possible. This treament is also useful to deal with the 
singularity of the global string in the next section. 

As inflation proceeds at the core of the string, the metric terms B, C, and H grow rapidly 
and the de Sitter expansion is established around the center of the string 

B C H /8ttG ta/ , " 

In addition to radial inflation, the string also inflates along the axial direction at the core. 
As shown in Fig. |5|, at the core, the ratio jr/% remains close to a constant (« 1). This 
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indicates that the expanding behavior along the z— direction is similar to that along the r— 
direction. 



III. GLOBAL STRING 

In this section, we consider the model of a U(l) global string. The action, field equations, 
and the metric are the same as those of the gauge string in the previous section without 
the gauge field (A^ = 0). In the absence of the gauge field the scalar field configuration of 
a global string stretches further radially. The situation is equivalent to that of the gauge 
string with (3 — > oo (e — > 0). On the basis of the cause of inflation discussed in the previous 
section, this guarantees that the global string will inflate at a lower symmetry-breaking scale 
and that the inflationary picture should not be much different. For n = 1, the critical value 
of 1] for inflation is found to be r) c ~ 0.23m p which is slightly lower than that of the gauge 
string. The critical value decreases as the winding number increases. The critical values are 
given for several winding numbers in Table |TJ. 

As mentioned in the Introduction, a static global string has a physical singularity outside 
the core. In this region, the energy-momentum tensor is given by 

T° = Tl = Ti = -T* = ^, (3.1) 

922 

and the closed-form solution of the Einstein equations was found by Cohen and Kaplan || 
ds 2 = (-) (~dt 2 + dz 2 ) + 7 2 (-) V2 (du 2 + d0 2 ) , (3.2) 



where uq = l/8nGr] 2 , and 7 is a constant of integration. This metric is transformed to the 



linearized solution obtained by Harari and Sikivie [ 10 1 when 77 m p 



ds 2 = [1 - 4G>(r)] (-dt 2 + dz 2 ) + dr 2 + [1 - 8G/i(r)] r 2 d6 2 . (3.3) 

Here /i(r) ~ JJ T^nrdr ~ 2irr] 2 ln (r/8) is the string mass per unit length out to a distance 
scale r from the core (~ 5). The relation between u and r is u ~ uq — In {r/5). Then the 



boundary of the core in the metric (|3.2j ) is located at u ~ uq and u = 00 corresponds to the 



center of the string. The solution ( |3.2|) exhibits coordinate singularities at u = and u = oo. 
The one at u = turns out to be a physical one and the existence of this singularity was also 
examined by analyzing the field equations in Refs. pT| , |T2|| . For Grand- Unified- Theory scale 
strings (rj ~ 10 16 GeV ~ lCU 3 m p ), the singularity is expelled outside the horizon, therefore, 
it does not cause any cosmological trouble. However, for supermassive strings, it is located 
near the core and could thus be problematic. 

In a similar way to the supermassive gauge string, one way to escape from this trouble 



is to introduce a time- dependent string. Indeed, Gregory |13| removed the singularity by 
taking a specific time- dependent form of the metric: only axial time dependence was intro- 
duced. But this metric still has a coordinate singularity (not a physical one). It was later 
questioned by Wang and Nogales |L7) that this coordinate singularity is unstable to both 
test particles and physical perturbations. In particular, they argued that the back reaction 
of the perturbations of null dust fluids will turn the coordinate singularity into a physical 
one. However, the situation is quite different for the supermassive strings. Because of the 
dynamical nature of the supermassive strings, we need to solve the system with a general 
time-dependent metric and scalar field. 

We now solve the equations numerically as we did for the gauge string. Fig. ^| shows 
the logarithmic values of the metric terms B, C, and H and Fig. |7| shows the Kretschmann 
scalar R a ^sR al316 at several moments of time. Their regular behaviors allow us to conclude 
that no singularity develops in the spacetime of a global string. For large i] global strings 
like supermassive gauge strings, the appropriate way to deal with their dynamical nature is 
to use a time-dependent treatment. 



IV. CONCLUSIONS 

We have investigated inflation in cosmic strings. In the core region, the strings undergo 
inflation radially as well as axially when rj > rj c . With unit winding number (n = 1) the 
critical values for inflation were found to be r] c « 0.25m p for a gauge string in the Bogomol'nyi 
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limit {{3 = 1) and r] c = 0.23m p for a global string. The critical values decrease as n and j3 
increase. We have explained this i] c variation in terms of the core size of defects. The core 
of defects inflates when its size becomes bigger than the horizon scale: for larger n and f3, 
strings have bigger cores, and the global string has a bigger core than the gauge string for 
a given 77. Regardless of the symmetry-breaking scale r], around the center of defects the 
de Sitter expansion is established since the scalar field stays about the top of the potenital 
(0 ~ 0). However, this is not sufficient for the cores of defects to inflate. Inflation requires 
another condition which is the core size being comparable to the horizon scale so that the 
core can be dynamical due to the gravitational effect. Or equivalently, the potential V(<p) 
needs to be flat enough at <fi ~ so that the field can spend enough time about the top of 
the potential. For this condition to be satisfied, the symmetry-breaking scale 77 needs to be 
sufficiently large. This description also explains why we have somewhat lower critical values 
of 7] for strings than those for domain walls and global monopoles (r] c « 0.33m p ). Strings 
have bigger cores at the same symmetry-breaking scale than the other defects. 

For supermassive gauge strings and all scale global strings, we have had troublesome 
physical singularities outside the core when we treat them in a static way. The elegant exit 
to nonsingular spacetimes is to introduce a time-dependent treatment. From the numerical 
simulations we could show that there is no singularity developing around time-dependent 
supermassive strings. 
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APPENDIX: FIELD EQUATIONS AND NUMERICAL ALGORITHMS 

Einstein's equations with the metric ( |2.3j ) and the energy-momentum tensor ( |2.5| ) are 
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G° = K\K\ + KlKl + K\K\ 

1 f C" H" B'C C'H' 

~~c IT ' 77 T 7 ~ TIT 



b 2 \ 

= 8ttG 



TT~b 
1 



12 , ^ , n2 Jlfl \ 2 . ^ 2 





2^ 


--) 


fi7 ~ 


CV 


Hr J 


(d 2 + 


a' 2 \ 
B 2 j 


+ v(<p) 



(Al) 



\(-G\ + G\ + G\) = K\ - (Kir ~ \{K\Kl - K\K\ + K\K\) 



+ 



2B 2 
AtiG 



'C" H" 

T + JT 



k'2 



b[_o_ 
~bT 

n 2 



a h' 
TIT 



3 B 2 + C 2 r 2 



<t> 2 (l-a) 2 - 



~TT~b 

n 2 1 
2e 2 C 2 r 2 



b' a h' \ 

h 2 

Br Cr Hr j 

/ « 2 + ^)-W) 



(A2) 



G\ 



+ 



(K 2 2 ) 2 - \(K\Kl + K 2 2 Kl - K\K\) 

'C" H" B'C' C'H' H'B' B' C' H' 
2B 2 \T~1T~ ~BT + TIT + ~HB~ ~ Br~ + 2 Tr~ + Hr , 



K 2 
1 



= 4nG 



12 



X 2 _l & 

+ W 2 



\2 n 

a) 2 + 



2e 2 C 2 r 2 



3a 2 



a 



/2 N 



(A3) 



2 {Gl + G 



+ 



2B 2 
AtiG 



(Kl) 2 - L -{-K\K 2 2 + K 2 2 Kl + KlK{) 

( C" H" B'C' C'H' H'B' B' C H' 
[~T + IT + ITT + TIT ~ JT~B + Br~~ ~Ch- + TTr~, 



k'2 



n 



i 



« 2 - 3^ I - V(cf>) 



a 



2e 2 C 2 r 2 



(A4) 



Goi = K 2 2 ' + Kf - (Kl - K 2 2 ) + ;) " (*l " *S)f 
n 2 1 

= ( A5 ) 



where 



^ = -|. ^ = -§' ^ 3 = -f- (A6) 
The field equation for <f> is 

+ ^ + + -^-|-^' + ^(i-«) 2 + > 2 -. 2 ) = o. 

(A7) 
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The field equation for a is 

a-(Kl-Kl + K!)a-^ + ^ 2 (^ + ^-^ + iy' + 2e^(a-l)=0. (A8) 

The field equations for the global string are obtained by setting a = 0. At t = 0, we assume 
a flat spacetime, 5(0, r) = C(0,r) = H(0, r) = 1, and zero velocities of scalar and gauge 
fields, 0(0, r) = d(0,r) = 0. Then we solve the scalar and vector field equations ( |A7| ) and 
(|A8|) to obtain 0(0, r) and a(0,r). K\{Q,r) and K$(0, r) are evaluated by the Hamiltonian 



and momentum constraint equations (|A1|) and (|A5[) after setting itf(0, r) = 0. In the 
next time step, B(t,r), C(t,r), and H(t,r) are calculated by Eq. (|A6|) . The Kl(t, r)'s are 
calculated by Eqs. ( [A2|) -(|A4]), and 0(t, r) and a(t,r) are calculated by Eqs. (|A7|) and ( |A8|) . 
The boundary conditions for and a are given in Eq. ( |2.4|) . We also use the boundary 
condition H'(t,0) = (Iff' = 0) and apply regularity condition Kl(t,0) = i^|(t, 0) at 
the origin. To solve the field equations numerically, we use the modified Crank-Nicholson 
method. We insist that the Hamiltonian and momentum constraint equations, ( |A1| ) and 
), are satisfied at every moment although they are not solved directly. 



We have only two constraint equations to set initial three if*(0,r)'s. So we assume 
fff(0,r) = which is consistent with the Einstein's equations. Nonetheless, it might not be 
the best choice. However, even if this choice is not the optimal one to describe our physical 
system, the numerical system will find its correct route propagating with the speed of light 
from the center of the string. The region in which we are interested relaxes to the correct 
configuration during the period of our numerical iteration. 
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TABLES 



n 


1 


2 


3 


4 


Vc/mp 


0.25 ± 0.005 


0.17 ±0.005 


0.13 ±0.005 


0.11 ±0.005 


(a) 


P 


0.5 


1 2 


3 


4 


Vc/rrip 


0.255 ±0.0025 


0.25 ±0.0025 0.24 ±0.0025 


0.235 ± 0.0025 


0.235 ± 0.0025 



(b) 

TABLE I. (a) The critical values of r\ for the gauge string in the Bogomol'nyi limit (3=1. (b) 
The critical values of rj for n = 1. 



n 


1 


2 


3 


4 


5 


Vc/mp 


0.23 ± 0.0025 


0.155 ±0.0025 


0.115 ±0.0025 


0.095 ± 0.0025 


0.075 ± 0.0025 



TABLE II. The critical values of rj for the global string 
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FIGURES 

FIG. 1. A plot of S/Hq 1 vs. n for the rj = 0.1, 0.5, 1.0, 2.0, 3. 0m p gauge strings (from the 
bottom to the top) in the fiat spacetime. The dashed line corresponds to sjn as a reference. Hq 1 
is the horizon size at the center of the string. 

FIG. 2. The scalar field configurations as functions of the proper radius CHqt at Hot = 0, 2, 4 
(from the left to the right) for the r\ = 0.5m p gauge string (n = 1, (3 = 1). The rapid growth of the 
proper radius in the core region indicates inflation. 

FIG. 3. Plots of (a) log 10 B, log 10 C, and (b) log 10 H vs. Hor at Hot = 1, 3, 5 (from the bottom 
to the top) for the rj = 0.5m p gauge string (n = I, (3 = 1). The metric terms behave regularly in r 
and t. 

FIG. 4. A plot of R a/3lS R a ^ s /H^ vs. H r at H Q t = 2, 4, 6 (from the right to the left) for 
the rj = 0.5m p gauge string (n = 1, f3 = 1). The scalar invariant is finite everywhere. The generic 
picture is not very different for different r/'s. 

FIG. 5. A plot of vs. H r at H t = 2,3,4,5 (from the bottom to the top) for the 

r\ = 0.5m p gauge string (n = 1, [3 = 1). The ratio stays close to a constant (« 1) in the core region 
(H r < 0.5). 

FIG. 6. Plots of (a) \og 10 B, log 10 C, and (b) log 10 H vs. H$r at Hot = 1,3,5 (from the 
bottom to the top) for the r] = 0Am p global string (n = 1). The metric terms behave regularly in 
r and t. 

FIG. 7. A plot of R aMS R a ^ 5 / 'H$ vs. H r at H t = 2, 4, 6 (from the right to the left) for the 
rj = 0.4m p global string (n = 1). The scalar invariant is finite everywhere. The generic picture is 
not very different for different 77's. The core is located at Hor < 0.5. 
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